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Abstract. The spin degree of freedom can play an essential role in determining the 
electrical transport properties of spin-polarized electron systems in metals or semicon- 
ductors. In this article, I address the dependence of spin-subsystem chemical potentials 
on accumulated spin-densities. I discuss both approaches which can be used to measure 
this fundamental thermodynamic quantity and the microscopic physics which deter- 
mines its value in several different systems. 



1 Introduction 

The role of the electronic spin degree of freedom in theories of the electrical 
transport properties of paramagnetic metals is passive; usually it as appears 
only as an afterthought — a factor of two to account for spin degeneracy. All this 
changes profoundly in electronic systems with substantial spin-polarization, ei- 
ther spontaneous or induced, particularly so if the system is either electrically 
or magnetically inhomogeneous. Recently, interest in the role of the electronic 
spin has increased, in part because of the possibility of fabricating technologi- 
cally useful magnetoresistive sensors and other devices based on spin-dependent 
transport effects, particularly giant magnetoresistance [1] and tunnel magne- 
toresistance [2] . Spintronics [3] , the study of spin-dependent electronic transport 
effects in systems containing metallic ferromagnets, is now a large and active 
area of basic and applied physics. In this article, I discuss the possibility of using 
transport experiments not to make devices, but instead to measure a fundamen- 
tal thermodynamic property of a spin-polarized electron system, the dependence 
of spin subsystem chemical potentials on accumulated spin-densities. This quan- 
tity can be important in modeling some spin-dependent transport effects. I will 
discuss several examples where it also provides a new and useful test of our 
understanding of the microscopic physics of a spin-polarized itinerant electron 
system. 

In Section 2, 1 derive a formally exact expression for the dependence of chemi- 
cal potentials on subsystem densities by defining a spin-dependent thermodynamic- 
density-of-states matrix. For bulk three-dimensional charged particle systems, 
the total density is fixed by electroneutrality requirements. In this case only the 
dependence of up and down spin chemical potentials on the difference between 
up and down spin densities is of interest. I show that these two quantities can be 
expressed in terms of the differential magnetic susceptibility, and a less familiar 
quantity, the derivative of chemical potential with respect to external field. The 
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focus of the article is this latter quantity, which I will refer to as the inverse 
magnetic compressibility . In Section 3, I discuss two spintronics experiments 
which can be exploited to measure its value in particular systems. In Section 4, 
I discuss the microscopic physics which determines its value in three distinct 
spin-polarized electron systems. Section 5 contains a brief summary. 



2 Thermodynamic Density-of-States Matrix 

Non-equilibrium spin accumulation [4, 5] due to electronic transport occurs 
generically in inhomogeneous spin-polarized electron systems and is a ubiqui- 
tous feature of spintronics. Any theory of spin accumulation requires, explicitly 
or implicitly, a model for the relationship between the up and down spin densi- 
ties and their chemical potentials. Linearizing around the equilibrium state, we 
can write 

d n = F U d n] + (DJ 1 + F U ) dn t . (1) 

where \i a is the spin-cr chemical potential and n a is the density of spin-cr elec- 
trons. In these equations, I treat nj and as separate thermodynamic vari- 
ables, something which is useful in discussing spin-accumulation since the pro- 
cesses which establish equilibrium between spin-f and spin-J. subsystems are 
often slow. The spin-quantization axis has been chosen to lie along the direction 
of net spin-polarization [6]. 

It is normally convenient to measure the local chemical potential of a charged 
particle system from the local electrostatic potential and this common convention 
is implicit throughout these notes. Accordingly 

_ 1 dF(T,n h n L ) 

with the free-energy per volume, F/V, calculated excluding any electrostatic 
contributions. The matrix of coefficients in Eq. (1) is accordingly given by the 
matrix of second derivatives of F/V with respect to n„, with the electrostatic 
term (which would diverge because of the long-range of the electron-electron 
interaction) neglected. (Note that F^[ = i*x,j.) To make contact with famil- 
iar descriptions of spin accumulation, I have introduced the band-theory spin- 
dependent density-of-states per volume, D a . If contributions due to correlation 
effects were neglected, the relationship between chemical potentials and densi- 
ties would be diagonal in spin indices and only the density-of-states terms would 
appear on the right-hand-side of Eq. (1). The density-of-states contribution to 
the spin-cr chemical potential change is simply the result of changing the filling 
of the spin-cr energy band. I will refer to the additional correlation terms in 
Eq. (1), F a ^>, as local-field corrections. The fact that the local-field corrections 
are in general off-diagonal in the spin indices can have a qualitative importance. 



Spintronic Spin Accumulation and Thermodynamics 



3 



Eq. (1) is formally exact, however, the values for the local field corrections, are 
not known in general, and their computation is a challenge to theory. In the 
following Section, I discuss two experiments which measure a particular combi- 
nation of these coefficients. 

Eq. (1) can be inverted to express spin-dependent density changes in terms 
of spin-dependent chemical potential changes 

dn„ = y^P cr<0 ./rf/x .' (3) 

where 

„ _ D a {l + D 9 F 9t9 ) 



1 + D a F a , a + D S F S . S - D a D s F^ s 

1 + D a F a>a + D s F StS - D a D s Fl s ' V ' 

where a —I if a =| and vice-versa. V = J^a a ' ^a,<r' ; the ra te of change of total 
density with chemical potential when the spin-subsystems are in equilibrium, is 
the total thermodynamic density-of-states of an electron system. If interactions 
are neglected V — + . In spintronics it is useful to generalize this concept 
by defining a thermodynamic density-of-states matrix as in Eq. (4). 

To evaluate the inverse magnetic compressibility, I add to the Hamiltonian 
a Zeeman coupling [7] term which contributes —g* fj, B H(n^ — njJ/2 to the free 
energy per unit volume. Here g* is the system's g- factor, fi B is the electron 
Bohr magneton, and H is the field strength. It is convenient to use a notation 
where \i„ is defined as the chemical potential without its Zeeman contribution, 
whereas [i is the full field-dependent chemical potential. Then the condition for 
equilibrium between up and down spins is 

M = Mt - 9*^bH/2 = m + g*fi B H/2, (5) 

or differentiating with respect to field strength: 

dn <9mt 1 _ dm + 1 



dg*fi B H dg*fi B H 2 dg*fi B H 2 

At fixed total electron density (dn-\ + dn^ = 0), differentiating Eq. (3) with 
respect to g* \i B H gives 



d^/dH _ _V M + V M _ D^+F 1A -F hl 
dm/dH~ P TiT +P iiT I)- 1 + F u - F u ' 

Combining Eq. (7) and Eq. (6), I obtain that for fixed total density 

cV £>j 1 - D~i 1 + Ffj - F^i 



dg*n B H 2[D- 1 + F ]A +F Li -2F 



(7) 



(8) 



UJ 
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Similar considerations lead to the following expression for the differential mag- 
netic susceptibility, 

= ffVg d(n T - n x ) = (ff>i?) 2 (Q] 

XS 2 dH D- 1 +Dl 1 +F 1 ,i+F u -2F hl - {) 

Both of these expressions are formally exact. 1 show in the following paragraph 
that the dependence of the chemical potentials on spin accumulation is specified 
by the inverse magnetic compressibility and xs- 

Theories of spin accumulation contain in general three elements: i) a the- 
ory of the spin and space dependent transport coefficients which lead to non- 
equilibrium spin-densities, ii) a theory for the disequilibration (/^ — fJ-i) pro- 
duced by these spin densities, and iii) a theory of the relaxation process which 
attempts to establish equilibrium between the spin subsystems. The thermody- 
namic property we are discussing is related to the second element. Assuming 
electroneutrality (dni = —dn^ ) it follows from Eq. (1) that 

_ -Pf 1 + F hl ~ F hl 
d(n-[ —ni) ^ 

dpi _ Dj 1 + F l,l - F t,l 
d ( n T- n l) 2 



(10) 



Note that d(^ — ^i)/d(n^ — nj.) = (g* Ms) 2 I^Xs- This relationship should 
not be a surprise since an external magnetic field maintains a chemical po- 
tential difference — Ml = 9* PbH and induces a magnetization per volume 
m = xsH = g*fi B (n^ — n{)/2. The linear relationship between /ij — ^ and the 
non-equilibrium spin accumulations is thus completely characterized by xs- It is 
only if we want to know the chemical potential shifts of up-spin and down-spin 
subsystems individually that the inverse magnetic compressibility is required. 
The individual chemical potential shifts driven by a non-equilibrium spin accu- 
mulation are: 



{g*H B ) 2 d( n] - ni ) d{g*fi B H) D" 1 + DJ 1 + F TiT + F u - 2F U 

4xs dm = _ 1 + 2 _ dn DJ 1 + F i,i — F hl 



(g*l^ B ) 2 d(n^ -n;) d(g*fi B H) D' 1 + DT 1 + F TiT + F u - 2F 



3 Spintronic Thermodynamic Measurements 
3.1 Field-dependent Coulomb Blockade Peaks 

The first type of experiment I discuss was pioneered by Ono and co-workers 
[8] and takes advantage of the equally spaced conductance peaks which occur 
in Coulomb blockade devices. The experimental geometry is that of a single- 
clcctron-transistor (SET) in which current flow from source to drain through 
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a small metallic particle is influenced by a gate voltage. In general source lead, 
drain lead, and metallic particle can be either paramagnetic or ferromagnetic. For 
dcfinitcncss, I assume that the only the small metallic particle is ferromagnetic 
and address a situation which is simpler than what has been encountered in 
experiments by also assuming a single domain. The following paragraph provides 
a simplified explanation the operation of a SET which is sufficient for present 
purposes. 

The dependence of the ground-state energy of an isolated metallic grain on 
its net charge is dominated by an electrostatic contribution and has the form 

e 2 (N - Nn) 2 

Eo(N) = { 2C °> + Ve (N/V) (11) 

where e(N — Nq) is the net charge on the grain and the effective capacitance of 
the grain C <~ R where R is the grain diameter. Here eo(n) is the energy per 
unit volume V calculated for a macroscopic grain which is electrically neutral. 
The conductance for current flow through the grain between lead electrodes is 
sharply peaked when the addition energy of the island is equal to the chemical 
potential in the lead electrodes i-e. when 

H(N) = E (N + 1) - E (N) = Mi . (12) 

In a SET the chemical potential for electrons on the metallic grain is manipulated 
by a gate voltage U: 

Eo(N) -> ^^^^ + Ve (N/V) - NeU. (13) 

As U is varied, the equilibrium number of particles on the grain changes. The 
number of particles in the grain's ground state changes between N and N + 1 
when Eq.( 12) is satisfied: 

eU^ = e 2 (N + l/2)/C + n(N/V)-fi L . (14) 

A peak occurs in the source-drain conductance at this value of the gate voltage. 
In this equation n(n) is the bulk chemical potential of an electrically neutral 
system, the quantity examined in Section 2. 

As discussed in Section 2, fj,(n) is field-dependent in general. It follows from 
Eq. (14) that provided the field-dependence of the lead chemical potentials can 
be ignored 

Mn) = d(eU* N ) 
d(g*fi B H) d{g*HBHy 1 ' 

the quantity of interest can simply be read off the gate voltage dependence of 
the Coulomb blockade conductance peak. More realistic models are complicated 
by geometry dependent cross capacitances between different circuit elements 
which invalidate the simple relationship between the chemical potential on the 
ferromagnetic grain and the gate voltage assumed here. However these two quan- 
tities are still proportional and the proportionality constant can be sorted out 
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experimentally by measuring the spacing between Coulomb blockade conduc- 
tance peaks at fixed voltage. Estimates of d{na{n) / d(g* hbH) have already been 
obtained [8] using this approach for one ferromagnetic transition metal. 

3.2 Field-dependent double-layer compressibility measurements 

The second potential experiment I discuss is a variant on one which has been 
used in the past [9, 10, 11, 12] to measure the compressibility of two-dimensional 
electron gas layers. It exploits techniques which have been developed [13] to make 
separate contact to nearby two-dimensional electron layers. The experimental set 
up can be thought of as a parallel plate capacitor, where one plate is a metal 
or heavily doped semiconductor and the second plate consists of two separately 
contacted two-dimensional electron layers, one on top of the other and closer to 
the metallic gate. A change in the charge density on the surface of the metal 
induces an opposing charge density distributed between the two two-dimensional 
electron layers. The equilibrium condition which determines the distribution of 
charge is: 

fj, T (n T ) = VB(n B ) + 4TTe 2 d(n - n T - n /2). (16) 

Here tit is the areal density in the top two-dimensional layer, ub is the density 
in the bottom two-dimensional layer, d is the separation between the two layers, 
and n = tit + riB is the total density in the two-layers. The second term on the 
right hand side of Eq. (16) is the electrostatic potential drop due to the electric 
field which exists between top and bottom two-dimensional layers. Note that the 
electric field between the gate and the top layer has magnitude Aire(n — no/2), 
i.e. en equals the surface charge density on the metallic gate up to a constant. If 
the top two-dimensional layer is held at ground, the gate voltage Vq is therefore 
proportional to n. The magnitude of the electric field below the bottom two- 
dimensional layer, 47ren /2, does not change during the experiment. 

To determine the compressibility [9], it is necessary only to measure the 
current which flows to the bottom layer when the gate voltage Vq (or equivalently 
n) is changed. Differentiating Eq. (16) with respect to n with ub + nx = n 
leads to [9] a relationship between a fundamental thermodynamic quantities 
and a conveniently measurable dimensionless experimental quantity exploited in 
previous experimental work: 

dn B _ du T /dn T 

dn 4:Tre 2 d + d[ir/dnT + dfis/dns 

Since the first term in the denominator of Eq. (17) is normally dominant the ex- 
periment provides a direct measurement of the dependence of chemical potential 
on density, inversely proportional to the compressibility, for the top layer. 

I propose using the same experimental setup to measure the dependence of 
the charge in the bottom layer on magnetic field at fixed gate voltage. Such a 
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measurement can be related to the inverse magnetic compressibility; differenti- 
ating Eq. (16) with respect to field at fixed n I find that 

dn B _ did T /d(g*fi B H) - d/j, B /d(g*fi B H) 
d(g*fi B H) 4-rre 2 d + dfj-r/driT + dfi B /dn B 

Since the denominator is dominated by the first term, the experiment provides 
a direct measurement of the difference between the rate of change chemical 
potential with field in top and bottom 2D layers. 



4 Microscopic Theory of the Thermodynamic 
Density-of-States Matrix 

4.1 Band Ferromagnets 

In order to contextualize the issues raised by these relatively new measure- 
ment possibilities, I first discuss the estimate for d/j,/d(g* /j, b H) which follows 
from the simplest possible mean-field theory of a ferromagnetic metal, Stoner- 
Wohlfarth [14] theory. In modern work this approach is wrapped in the cloak of 
spin-density-functional theory [15]. In the Stoner-Wohlfarth theory, quasiparti- 
cle spin-up and spin-down energies for each band and for each wavevector in the 
crystal's Brillouin-zone are split by an amount 

A = 2/j, B H eff = 2(Im + fi B H). (19) 

The effective magnetic field includes an exchange contribution which is propor- 
tional to the magnetic moment per Bohr magneton per volume of the system, 

m= =nt-rii. 20 

VfJ, B 14 

where M = fj, B (N^ — N±) is the total magnetic moment. The exchange integral 
J is a phenomcnological material property which is characteristic of a given 
system. In a given field, m and the chemical potential \i are determined by self- 
consistently occupying spin-split bands. At T = the self-consistent mean-field 
equations are 

fH+A/2 

m= I deD(e) 

Jfi-A/2 

pH-A/2 rfi+A/2 

n = 2 deD(e)+ / deD(e) (21) 

J-oo J ix— A/2 

where D(e) is the density-of-states per volume per spin at A = 0. Let A and I 
denote the self-consistent exchange splitting of the band and the self-consistent 
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value of m in the absence of an external field respectively. Expanding \i and A 
to first order in if, I find that 

da D7 1 - D7 1 

^ (22) 



d(g*u B H) 2[D~ i l + D- 1 - AI] 

Here = D(uq + Ao/2) and D\ = D(/xo — Aq/2) are the majority and minority 
spin densities of states in the absence of a field. 

This expression can be recognized as a special case of the general result 
derived in Section 2 which is obtained by holding the bands rigid and taking 
F^i — — = — -Fj.4 = /. In the formally exact theory, the three local-field- 
factors are independent. Density-functional theory [15] is a practical approach to 
many-particle physics which has been applied successfully to evaluate the ther- 
modynamic properties of metals and underpins the modern theory [16] of ferro- 
magnetic transition metals. In density-functional theory, the local field factors 
arise from an interplay between band-structure details and exchange-correlation 
single-particle potentials. Expectations based on the electron gas case, discussed 
for two-dimensions below, suggest that the diagonal local-field-factors should 
be negative and larger in magnitude than the positive off-diagonal local-field- 
factor. It is not at all obvious that these expectations apply to transition met- 
als, especially because the majority spins at the Fermi energy tend to have 
predominantly itinerant s-electron character while the minority spins tend to 
have predominantly localized cZ-electron character. Comparison of theoretical 
and experimental values for both the differential magnetic susceptibility, and 
the magnetic compressibility, presents a serious and interesting challenge to the 
density-functional theory of metallic magnetism. 



4.2 Zero Field Two-Dimensional Electron Gas 

The ground state of the two-dimensional electron gas is not ferromagnetic except 
at extremely low densities [18, 19]. Spin-polarization can, however, be induced by 
application of an external magnetic field [20]. I now discuss a simple theory for 
the value oidfi/d{g* ^bH) in a spin-polarized two-dimensional electron gas which 
is based on the Hartree-Fock approximation [21] for its energy: two-dimensional 
electron system: 

E nuni nj+nj 8e>f + nf] 

A 2D 2D 37^/2 1 > 

Here D a = m* /(2nh 2 ) is the density-of-states per area A in a 2D system. The 
first term on the right hand side of Eq.( 23) is the band energy and the second 
term is the exchange energy. Differentiating twice with respect to density we 
obtain for the diagonal local-field-correction 

^ - (24) 
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The off diagonal local-field-corrections vanishes in this approximation because of 
the neglect, in the Hartree-Fock approximation, of correlations between electrons 
of opposite spin. In the presence of an external magnetic field the partitioning 
of density between spin-subsystems is determined by setting fj,^ — fj,^ = g* [i B H 

where \x a = n a / D2D — ^e 2 n]J 2 /i: 1 / 2 . For electron-gas density[17] parameter r s = 
2, the Hartree-Fock ground state is paramagnetic and the spin-polarization is 
initially linear in field. At stronger fields, the exchange term causes the minority- 
spin chemical potential to increase with decreasing density. The polarization 
must then increases more rapidly with field in order to establish the required 
chemical potential difference. The polarization eventually jumps from a partial 
value to full polarization when g* /.ibH ~ 0.034Ry. (The atomic energy unit 
lRy = e i m*/2e 2 h 2 is w 5.5meV in GaAs.) 




Fig. 1. Hartree-Fock approximation majority-spin (solid line) and minority-spin 
(dashed line) chemical potentials vs. Zeeman coupling strength for a two-dimensional 
electron gas with density parameter r s = 2. 



The results for the Zeeman coupling denpendence of subsystem chemical po- 
tentials and spin-polarization, shown in Fig.( 1) and Fig.( 2) respectively, can 
lead to substantial values of dfi/d(g* hbH) which can even exceeds one over a 
wide range of Zeeman coupling strength. This is despite the fact that the minor- 
ity and majority spin density-of-states are identical in two-dimensions so that the 
inverse magnetic compressibility would be identically zero for a non-interacting 
electron system. Numerical results for the inverse magnetic compressibility for 
this illustrative example arc shown in Fig.( 3). Of course, we do not expect the 
Hartree-Fock approximation theory to be reliable in an electron gas for r s > 1. 
The large spin-polarizations at relatively weak Zeeman couplings in Fig. (2), oc- 
cur because the Hartree-Fock approximation ground state is ferromagnetic for 
r s > 7I-/2 1 / 2 . We know from reliable quantum Monte Carlo calculations [18, 19] 
that, in reality, this instability does not occur until much lower densities are 
reached. Nevertheless illustrative calculation does suggest what should be ex- 
pected in low-density two-dimensional electron gas layers. For modern samples, 
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g u B H [Atomic Units] 

Fig. 2. Hartree-Fock approximation m/n= {n-\ — rti)/{n-\ + n{) for a two-dimensional 
electron gas with density parameter r a — 2 vs. Zeeman coupling strength g* fisH in 
atomic units (e 4 m* /2h 2 ). For GaAs the Zeeman coupling strength in atomic (Ry) units 
is 0.0052// [Tesla]. The system jumps to full polarization at g* i±bH ~ 0.034Ry. 
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g u B H [Atomic Units] 

Fig. 3. Hartree-Fock approximation d/j,/d(g* hbH) vs. Zeeman coupling strength, for 
a two-dimensional electron gas with density parameter r s = 2. For g* hbH > 0.034, 
only the majority spins are occupied and d/j,/d(g* HbH) = —0.5. 

which have high mobility at densities a few xl0 lo cm~ 2 , it is clear that a sub- 
stantial degree of spin-polarization (and possibly complete spin polarization) can 
be achieved in laboratory superconducting magnets. Measurement of the inverse 
magnetic compressibility, coupled with simultaneous measurements of the spin 
polarization, will provide a strict test of our understanding of magnetism in 
strongly correlated electron gas systems. 

4.3 Quantum Hall Ferromagnets 

In the quantum Hall regime, the 2D electron density is measured in units of the 
density which can be accommodated by a full Landau level; the Landau level fill- 
ing factor v = (27r£ 2 )n where £ = (fic/eB) 1 / 2 is the magnetic length. At Landau 
level filling factor v = 1 , the 2D electron system has a strong ferromagnet ground 
state, i.e. all spins can be aligned with the field by arbitrarily weak Zeeman cou- 
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pling. Quantum Hall ferromagnets have unusual charged excitations known as 
Skyrmions [22, 23]. Skrymions are topologically non-trivial configurations of the 
magnetization orientation distribution which occur in any 2D magnet, but carry 
an electrical charge e only in the case of quantum Hall ferromagnets. Skyrmions 
have an internal [24] integer quantum number K which specifies the number 
reversed spins in their interior. When the Landau level filling factor is close to 
v = 1 , the low energy states of the system can be described in terms of Skyrmion 
degrees of freedom. 

A simple model of a 2DES, valid at low temperature in this regime because 
Skyrmions are dilute, is obtained by ignoring Skyrmion-Skyrmion interactions. 
The following grand-canonical ensemble expressions specify the occupation prob- 
abilities of the Nff, = A/(2tt£ 2 ) Skyrmion quasielectron and quasihole states with 
K excess reversed spins: 



Here /(e) = (exp(e/fcsT) + is a Fermi factor [25], €k is the energy of a 
Skyrmion quasiparticle, (2tt£ 2 )~ 1 is the density of a full Landau level, and we 
have chosen the zero of energy so that quasielectron and quasihole skyrmion 
states have the same energy [26]. When the spin subsystems are in equilibrium, 
we can use Eqs. (25) to calculate the chemical potential, given the Landau level 
filling factor. The Landau level filling factor is increased by quasielectron exci- 
tations and decreased by quasihole excitations: 



Eq.( 25) follows from the property that formation of the K — th quasielectron 
Skyrmion requires the addition of K + 1 spin-down electrons and the removal 
of K spin-up electrons from the v = 1 ground state, while formation of the 
K—th quasihole Skyrmion requires the addition of K spin-down electrons and the 
removal of K + 1 spin- up electrons. For non- interacting electrons only the K = 
quasiparticles occur; for typical 2DES's, on the other hand, the lowest energy 
quasiparticlcs have K — 3; these quasiparticles dominate the low-temperature 
properties of the system for v close to 1. From Eqs.( 25) I obtain the following 
thermodynamic density-of-states matrix:, 



n Ke = f(e K + -(K+ 1)^) 
n K h = f{tK + (K + l)/x T - K/j,^. 



(25) 




(26) 



K 



P T , T = (27re 2 )- 1 ]T [(K + l) 2 A(e K + /i) + K 2 A(e K - M )] 



(27) 



K 



V hl = -(27T^ 2 )- 1 ]T K(K + 1) [A(e K + /x) + e K - //)] 



(28) 



K 



where A(x) = sech 2 (x/ '2) /4k B T. 

These expressions have been used previously to analyze [28] spin bottlenecks 
which have been observed [29] in the quantum Hall regime. Here I evaluate the 
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dependence of the equilibrium chemical potential on Zccman coupling strength. 
The possibility of measuring this quantity in two-dimensional electron systems 
was discussed in the previous section. Zeeman coupling adds [22, 26] to the 
Skyrmion quasiparticle energies, — > £k + 9* HbH{K + 1/2). In order for the 
filling factor to be held fixed as the Zeeman coupling strength varies, the chemical 
potential must change: 



c) M _ dn/d(g*^ B H)\^ _ j: K (K + l/2)[A(e K — fi) — A(e K + M )] 



d(g*fi B H) dn/dfi\ H J2 K [A(e K - n) + A(e K + rf} 

(29) 

We see in Fig. (4) that for an ideal disorder-free two-dimensional electron sys- 
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x 



k B T=0.02 



Fig. 4. dfi/d(g* hbH) vs. Landau level filling factor for a quantum Hall system near 
v = 1. at temperature ksT = 0.02e 2 /£. 



tern, d/j/d(g* [j,bH) is typically larger than one and changes sign when the filling 
factor crosses v = 1. These results were calculated using Skrymion state ener- 
gies evaluated by Palacios et. al. [27]. The value of this quantity gives infor- 
mation about the Skyrmion system complementary to that available from other 
experiments [23]. Its value depends both on the average K value of the Skrymions 
present in the system, the quantity revealed by Knight shift measurements, but 
also on the way in which the distribution of K values changes with field. Unlike 
transport experiments, this thermodynamics measurement would be sensitive to 
Skrymions localized by disorder in the system. 
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5 Summary 

In this article, I have discussed the dependence of chemical potential on exter- 
nal magnetic field in a spin-polarized electron system. This quantity appears 
directly in theories of electron systems disturbed by non-equilibrium spin accu- 
mulations generated by spin and space-dependent electronic transport. I have 
discussed techniques which can be used to measure this fundamental thermody- 
namic quantity, both in two-dimensional electron systems and in ferromagnetic 
metals. I believe that its measurement can provide an important and interesting 
test of theory in several quite different spin-polarized electron systems. 
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